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; ABSTRACT 
. 

^ ■ The magnetic field strength in molecular clouds is a fundamental quantity for 

^ ' theories of star formation. It is estimated by Zeeman splitting measurements in a 

■ few dense molecular cores, but its volume-averaged value within large molecular 

clouds (over several parsecs) is still uncertain. In this work we provide a new 
method to constrain the average magnetic field strength in molecular clouds. We 
compare the power spectrum of gas density of molecular clouds with that of two 
350^ numerical simulations of supersonic MHD turbulence. The numerical simu- 
lation with approximate equipartition of kinetic and magnetic energies (model A) 
yields the column density power spectrum P{k) oc /j;-2.25±o.oi^ ^^^e super-Alfvenic 
simulation (model B) P{k) oc /j;~2.7i±o.oi rjj-^g column density power spectrum of 
the Perseus, Taurus and Rosetta molecular cloud complexes is found to be well 
approximated by a power law, Po{k) oc k""-, with a = 2.74 ±0.07, 2.74 ±0.08 and 
2.76 ± 0.08 respectively. We conclude that the observations are consistent with 
the presence of super-Alfvenic turbulence in molecular clouds (model B) while 
model A is inconsistent (more than 99% confidence) with the observations. 
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1. Introduction 

The volume-averaged magnetic field strength in molecular clouds has never been mea- 
sured directly. Zeeman splitting has been detected only from a few dense molecular cloud 
cores, where emission lines from molecules such as OH and CN are observed (Crutcher 
1999; Bourke et al. 2001). Dense cores fill only a small fraction of the volume of molecular 
clouds. Therefore, the magnetic field strength averaged over the molecular cloud volume 
cannot be directly inferred from its value in the cores. This is true especially if the magnetic 
field strength has a very intermittent distribution and is correlated with the gas density, as 
suggested by Padoan & Nordlund (1999). 

Estimates of magnetic field strength in molecular clouds have been inferred from the 
dispersion in the polarization angle (e.g. Myers & Goodman, 1991; Chrysostomou et al., 

1994; Lai et al., 2001; Matthews & Wilson, 2002; Lai et al., 2002) as originally suggested 
by Davis (1951) and Chandrasekhar & Fermi (1953). This method was tested in numerical 
simulations of MHD turbulence by Ostriker et al. (2001), Padoan et al. (2001) and Heitsch 
et al. (2001). The relative motion of ions and neutral molecules, as manifested by a com- 
parison of their spectral lines, has also been used to estimate the magnetic field strength in 
molecular clouds (Houde et al. 2000, 2002). 

In this work we present a new way to constrain the average magnetic field strength in 
molecular clouds, based on the density power spectrum. In § 2 we show numerical simulations 
of supersonic MHD turbulence with different magnetic field strength yield different power 
spectra of gas density. The power spectrum is then computed from maps of molecular clouds 
in § 3. We find that only a rather low value for the average magnetic field strength, leading to 
super-Alfvenic turbulence, is consistent with the observations. Conclusions are summarized 
in § 4. 

The density power spectrum depends on the rms sonic Mach number of the turbulence, 
Ms (the ratio of the rms fiow velocity and the sound speed). In this work we use only 
simulations with Ms ~ 10 because that is the approximate value in the molecular cloud 
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complexes we have studied. We do not study the dependence of the power spectrum on the 
vahic of Afs. In general, smaller values of Ms yield steeper density power spectra than the 
Ms ~ 10 models (this was verified with a set of simulations that will be presented elsewhere). 
Turbulent flows with Ms ^ 1, for example, are expected to generate a Kolmogorov density 
power spectrum, proportional to k~^^^^. This slope is comparable with the electron density 
power spectrum on very small scale estimated from scintillation studies (Armstrong et al. 
1995). HI surveys of our galaxy and the Magellanic Clouds (Crovisier & Dickey 1983; 
Elmegreen et al. 2001; Stanimirovic & Lazarian 2001; Dickey et al. 2001) produced power 
spectra with slope intermediate between the present results in molecular clouds and the 
scintillation studies. 



2. Power Spectrum of Gas Density in Supersonic MHD Turbulence 

In order to study the power spectrum of gas density, we have run two simulations of 
driven supersonic MHD turbulence with rms sonic Mach number Ms ~ 10 and isothermal 
equation of state. We have solved the three dimensional compressible MHD equations in a 
staggered mesh with 350^ computational cells and periodic boundary conditions. The initial 
magnetic and density fields are uniform. The flow is driven by an external large scale random 
and solenoidal force, correlated at the largest scale turn-over time. The time derivative of 
the random force is generated in Fourier space, with power only in the range of wavenumbers 
1 < k Linesh/2 TT < 2 . The initial velocity field is proportional to the initial force, with an rms 
amplitude of approximately 50% of its relaxed value. Details about the numerical method 
are given in Padoan & Nordlund (1999). 

We have run the simulations for five dynamical times. The dynamical time is here 
defined as = Linesh/(2M), where u is the rms flow velocity. We choose this definition of 
dynamical time because the flow is forced up to wavenumber k — Air/L^esh and therefore 
the largest turbulent scale is approximately L^^sh/'^- 

We characterize the simulations based on the relative importance of magnetic and dy- 
namic (turbulent) pressure. The pressure ratio is defined as Pm/Pd = (-S^)/[87r (pw^)], av- 
eraging over the computational volume and over the last four dynamical times. In the first 
simulation (model A) the statistically relaxed pressure ratio is Pm/Pd ~ 0.65 ±0.05 (approx- 
imate equipartition of magnetic and kinetic energies). In the second simulation (model B), 
the ratio is Pm/Pd ~ 0.09 ± 0.01 (kinetic energy of the turbulence approximately five times 
larger than magnetic energy). The standard deviation of approximately 10% in the pressure 
ratio of models A and B is due primarily to time fiuctuations of the random force. 
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Power spectra have been computed for 18 times over the last four dynamical times 
(the small scale portions of these power spectra are statistically independent because the 
dynamical time decreases with spatial scale). The slopes are computed from a least square 
fit of the time-averaged power spectra plotted in Figure 4. Wc find that the power spectrum 
of the density field is sensitive to the pressure ratio (or the average magnetic field strength). 
The power spectrum of the three dimensional (3D) density field is P{k) oc /;;-2.25±o.oi ^j^g 
equipartition run (model A) and P{k) oc fc-2.7o±o.oi ^j^g super-Alfvenic case (model B). 

In isotropic turbulence the power spectrum of the projected density is the same as the 
power spectrum of the 3D density field (not necessarily in the presence of a mean magnetic 
field or in real molecular clouds). Wc have verified this in our numerical data. We have also 
verified that the power spectrum of the projected density does not depend on the direction 
of projection relative to the direction of the mean magnetic field. However, our results are 
based on the power spectra of the 3D density field and not of the projected density, because 
the statistical sample size is reduced by the projection (larger noise in 2D power spectra 
than in 3D ones). 



3. Power Spectrum of Column Density in Molecular Clouds 

The distribution of column density in molecular clouds can be estimated from maps of 
the J=l-0 ^^CO emission fine. In this work we use J=l-0 ^^CO maps of the Perseus (Padoan 
et al. 1999), Taurus (Mizuno et al. 1995) and Rosetta (Heyer & ct al. 2003) molecular 
cloud complexes and find the power spectrum of projected density in these regions is well 
approximated by a power law, Po{k) oc k~^°. 

The value of the gas column density inferred from the J=l-0 ^^CO maps depends on the 
distribution of ^^CO abundance and J=l-0 ^^CO excitation temperature, Tg^- The column 
density can be estimated using the LTE method (Dickman 1978; Harjunpaa & Mattila 1996; 
Padoan ct al. 1998a). In the LTE method the value of T^x along each line of sight is assumed 
to be constant and is estimated using the observed peak temperature, Tj. = Tp, of an optically 
thick line (r ^ 1) in the equation: 

Z^[J{T,^)-J{n,)]{l-e-^) (1) 

where Tbg = 2.7 K is the background temperature, and the function J{T) is defined as: 

with To — huio/k and is the frequency of the J=l-0 ^^CO transition. The J=l-0 ^^CO 
transition, when available, is generally used as the optically thick line to estimate the value 
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of Tex- For the present analysis we apply the LTE method using only the J=l-0 ^^CO line 
because ^^CO maps are not available to us. This line is optically thick only in the brightest 
regions of the maps used in the present work. The value of Tgx can be estimated from the 
J=l-0 ^'^CO line only in those regions. Therefore, we use the peak temperature of this line 
over the whole map to estimate the value of Tex as outlined above. This single value of Tex 
is used for all map positions (the effect of this assumption is addressed below with radiative 
transfer calculations). Finally, the ^^CO column density is given by: 

where t{v) is the optical depth in the velocity channel corresponding to the velocity v and is 
estimated from equation (1) using the estimated value of Tgx and the radiation temperature 
Tr(v) given by the observed line profile. Av is the width of the velocity channels in the 
observations, expressed in km/s, and the partition function Q is assumed to be a constant 
over the map so its value is not required as it does not affect the slope of the density power 
spectrum. 

Column density maps of the Perseus, Taurus and Rosetta regions have been computed 
from the observed J=l-0 ^^CO spectral maps with this simphfied LTE method. Their power 
spectra are plotted in Figure 4. They are shown to be well approximated by power laws 
over one decade or more in wavenumber, Po(^) k~^°, with Oq ~ 2.87 ± 0.04 for Perseus, 
Co ^ 2.87 ± 0.06 for Taurus and ^ 2.89 ± 0.06 for Rosetta. 

In order to estimate the effect of assuming constant Q and Tex, the simphfied LTE 
method for estimating the gas column density has been applied to synthetic spectra. As- 
suming a uniform ^^CO abundance, synthetic spectral maps of the J=l-0 ^^CO hue have 
been computed with a 3D non-LTE Monte Carlo radiative transfer code (Juvela 1997) as 
described in Padoan et al. (1998b). The radiative transfer is computed through the density 
and velocity field of two snapshots of model B (at 3 and 4 dynamical times), regriddcd to a 
resolution of 175^, assuming average density (n) = 500 cm^^ and mesh size Tmcsh = 10 pc. 
The 3D distribution of kinetic temperature is computed self-consistently as part of the ra- 
diative transfer solution, from the balance of cosmic ray heating and molecular and atomic 
cooling. We have assumed a cosmic ray hydrogen ionization rate of 2 x 10~^^ s~^ and heating 
of 8 eV per ionization. 

Six maps of 175 x 175 synthetic J=l-0 ^^CO lines are obtained in this way, corresponding 
to three orthogonal directions of projection of the two 3D snapshots. The simphfied LTE 
method is then applied to these six synthetic spectral maps and the estimated column density 
is compared with the actual column density in the original snapshots of model B. The result 
for one of the maps is plotted in Figure 4. Figure 4 shows the estimated column density. 
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A^LTE, is roughly proportional to the true column density, but tends to saturate at large 
column density values. This is due both to the low gas temperature and to the J=l-0 
^^CO line saturation in the densest regions. The power spectra of the column density field 
estimated from the synthetic spectral map, Plte(^) oc A;^'*^™, arc then compared with 
the power spectra computed directly from the projections of the original MHD data cube. 

We find that Plte(^) is generally steeper than Pmhd(^)- This is primarily due to 
the decreased emission from the densest regions due to their low gas temperature and to 
the saturation of the ,1=1-0 ^'^CO line. The difference is Clte — omhd = 0.13 ± 0.06. If 
this correction is applied to the observed molecular cloud complexes, the corrected power 
spectrum has a slope a = 2.74±0.07 for Perseus, a = 2.74±0.08 for Taurus and a = 2.76±0.08 
for . At the 1 a level, model B is consistent with these molecular cloud complex power spectra, 
while the power spectrum of model A is inconsistent with the observations with virtually 
100% confidence (7 a). 

Variations in ^^CO abundance may affect the estimated column density. However, sig- 
nificant CO depletion is expected only above 10 magnitudes of visual extinction and the CO 
abundance should drop only below 1-2 magnitudes. Most of the gas mass in molecular cloud 
complexes emits at values of visual extinction between 1 and 10 mag. Therefore, the effect 
of gas temperature variations estimated above should be more important than the neglected 
effect of variations in ^^CO abundance. 



4. Conclusions 

The average magnetic field strength in molecular clouds cannot be measured directly. 
However, it can be inferred from observational data due to its effect on the gas dynamics. 
In this work we have found the power spectrum of the density field is a sensitive diagnostic 
of the magnetic field strength. Numerical simulations of supersonic MHD turbulence with 
rms sonic Mach number Ms ^ 10 develop a power law power spectrum of gas density, 
-PmhdI^) oc A;~"mhd_ xhe value of the power law exponent is cmho = 2.25 ± 0.01 when 
the flow rms velocity is comparable to the Alfven velocity and aMHD = 2.71 ± 0.01 in the 
super- Alfvenic simulation. 

The density power spectrum can be measured also in molecular cloud complexes, for 
example using J=l-0 ^^CO maps. However, the column density (and its power spectrum) 
estimated using only the J=l-0 ^^CO line and the LTE method are biased by a number of 
uncertainties. The most signiflcant uncertainties are the 3D distribution of the gas kinetic 
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temperature in the molecular cloud complexes and the saturation of the J=l-0 ^^CO line 
in very dense regions. We have investigated the effect of these uncertainties in the density 
power spectrum using synthetic J=l-0 ^'^CO spectral maps computed with a non-LTE Monte 
Carlo radiative transfer code. The 3D equilibrium temperature distribution is computed 
self-consistently as part of the radiative transfer solution by balancing cosmic ray heating 
with molecular and atomic cooling. The correct power spectrum slope, aMHD, is found to be 
smaller than the slope estimated with the LTE method, oltE) with omho = olte — 0.13±0.06. 

With this correction, the power spectrum slope is a = 2.74 ± 0.07 for Perseus, a = 
2.74 ± 0.08 for Taurus and a = 2.76 ± 0.08 for Rosctta. The supcr-Alfvcnic model is 
consistent with this result, while the model with rms flow velocity comparable to the Alfven 
velocity is ruled out by the observations. This is yet another indication that super-Alfvenic 
turbulence provides a good description of molecular cloud dynamics. As first proposed by 
Padoan & Nordlund (1997, 1999), the average magnetic field strength in molecular clouds 
may be much smaller than required to support them against the gravitational collapse. 
Evidence of super-Alfvenic turbulence was also found recently by Troland and Heiles (2001) 
in HI clouds. 

MJ acknowledges the support from the Academy of Finland Grants no. 174854 and 
175068. Computing time was provided by the Danish Center for Scientific Computing. The 
research of RJ is partially supported by NSF grant AST-0206031. 
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Fig. 1. — Time-averaged density power spectra of model B (upper plot), and model A 
(lower plot). The solid lines are least square fits in the range 0.3-0.9 pc~^ for model A and 
0.3-1.2 pc~^ for model B, assuming a mesh size Lmesh = 10 pc. 
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Fig. 2. — Power spectra of three molecular cloud complexes. The power spectra are computed 
from images of column density obtained with the LTE method applied to maps of the J=l-0 
^^CO emission line. 
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Fig. 3. — Scatter plot of column density estimated with the LTE method applied to a 
synthetic map of the J=l-0 ^^CO emission line versus the column density in the MHD data 
cube used to compute the synthetic spectra. 



